Numerical Time Integration for Air Pollution Models

Due to the large number of chemical species and the three space dimensions, o -the-shelf sti ODE integrators are not feasible for the numerical time integration of sti systems of advection-diffusion-reaction equations @c @t +r (uc) = r (Kr c) +R (c) ; c = c(x; t); c 2 IRm; x 2 IR3 from the eld of air pollution modelling. This has led to the use of special time integration techniques. This paper is devoted to a survey of such techniques, encompassing sti chemistry solvers, positive advection schemes, time or operator splitting, implicit-explicit methods and approximate matrix factorization solutions. Of great importance in practice is high performance computing due to the huge problem scales, in particular for global models. We will therefore also report on experiences with vector/parallel shared memory and massively parallel distributed memory architectures and clusters of workstations. The survey is not entirely unique to air pollution models and biased towards work done at CWI over approximately the last 5 years. 1991 Mathematics Subject Classi cation: Primary: 65M06, 65M20. Secondary: 65Y05, 65Y20. 1991 Computing Reviews Classi cation System: G.1.8 and G.1.1, J.2.